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Abstract. We discuss a two-parameter family of maps that generalize piece- 
wise linear, expanding maps of the circle. One parameter measures the effect 
of a non-linearity which bends the branches of the linear map. The second 
parameter rotates points by a fixed angle. For small values of the nonlinearity 
parameter, we compute the invariant measure and show that it has a singular 
density to first order in the nonlinearity parameter. Its Fourier modes have 
forms similar to the Weierstrass function. We discuss the consequences of this 
singularity on the Lyapunov exponents and on the transport properties of the 
corresponding multibaker map. For larger non-linearities, the map becomes 
non-hyperbolic and exhibits a series of period-adding bifurcations. 

1. Introduction. The concept of natural invariant measure is central to the the- 
ory of Sinai-Ruelle-Bowen measures. Rigorous results pertaining to the existence 
and uniqueness of such measures have been obtained for the most part in the frame- 
work of Anosov diffeomorphisms and Axiom A systems. However, as discussed by 
Chernov £Q, there are hyperbolic maps other than diffeomorphisms for which SRB 
measures were constructed [181 1221 12]. It is our purpose to investigate some of the 
properties of a class of such systems with singularities. 

This paper is concerned with the parametric dependence of a particularly simple 
class of one-dimensional circle maps which may be regarded as non-linear gener- 
alizations of a class of piece-wise linear, expanding maps. The generalized maps 
have two essential properties, the first being that these maps are singular in the 
sense that their first derivative posseses a singularity (or a finite set of them), and 
second that the maps perform some irrational rotation such that the singularity 
is moved ergodically about the circle upon iteration. This class of maps will be 
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characterized by two parameters: a first parameter which measures the strength of 
the nonlinearity, and a second one which is a rotation angle. 

Depending on the strength of the non-linearity, the map can be hyperbolic or be 
non-hyperbolic. In the framework of the one-dimensional circle maps we consider 
here, these two regimes are best described by measures on interval-filling attrac- 
tors or by delta-measures, respectively. It is convenient to assume that the linear 
regime, i. e. when the non-linearity is set to zero, corresponds to a regime where 
the map is uniformly expanding, for instance, an angle-doubling automorphism. 
The onset of non-linearity perturbs this hyperbolic regime without altering this 
property. However, when the non-linearity reaches a threshold value, the map loses 
its hyperbolicity and stable periodic orbits may occur. 

The presence of a rotation has interesting consequences in both these regimes. 
In the regime of small non-linearity, the hyperbolicity will not be affected by the 
rotation since the map is everywhere expanding. However, because of the exis- 
tence of singularities, upon iteration, an initially uniform density will evolve into 
a discontinuous one. The invariant density has a dense set of discontinuities, but 
remains finite everywhere. For larger non-linearity, the rotation can take a non- 
expanding region into an expanding one, with the consequence that the attractor 
may alternate between hyperbolic and non-hyperbolic regions as the intensity of 
the non-linearity or the rotation angle are varied. 

2. Singular circle maps. Consider the family of automorphisms of the circle 

</Jo,c : x — > 2x + C mod 1 , (1) 

which combines a rotation by an angle £ with a uniformly expanding automorphism. 
It is clear that the Lebesgue measure is invariant under such a map, regardless of 
the parameter The parametric dependences of maps with a rotation become less 
trivial when the linear expansion is replaced by a nonlinear one. Consider maps of 
the form 

fax '■ x — > 2 fa{x) + C mod 1 , (2) 

where we assume f a is defined over the unit interval, with the following features, 
motivated by previous studies of non-linear baker maps^I] : 

(i) the end points 0, 1 of the interval are fixed points, / a (0) = and / a (l) = 1, 
such that the former is an attracting fixed point, f' a (0) < 1, while the latter 
is repelling, f' a {l) > 1; 

(ii) f a is analytic with respect to the parameter a; 

(iii) f a {x) is Holder continous with respect to the variable x. Moreover, all the 
derivatives of f a with respect to x exist at the end points; 

(iv) fa has a semi-group property, namely the composition of fa with itself yields 

/2a, 

fa ° fa{x) = f 2a {x) , (3) 

and the inverse of f a is obtained by changing the sign of a, 

f-\x) = f. a (x) ; (4) 

( v ) fa has the following property of reversibility 



/ (1 - z) = 1 - /_ a (z) 



(5) 
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The semi-group property, Eqs. II) . imply that f a is the identity for a = 0, 
namely f (x) = x; 

We note that Eqs. I)4I5[I . imply the following : 

mm = i ■ (6) 

Moreover these quantities are different from 1 as required in item (i). 

In this paper, we will consider in its generality the class of perturbations f a 
satisfying the above properties. In particular, we will consider the small parameter 
regime, \a\ <C 1, where f a can be written as a power expansion in its parameter, 

f a (x)=x + af (1 \x)+0(a 2 ), (7) 

and demand that the property of the derivative of /„ at the extremities of the 
interval, see item (i) above, hold to first order in a. 

An immediate property of follows from Eq. namely it is symmetric with 
respect to x — 1/2, 

fW(x) = fV>(l-x), (8) 

which also implies that the derivatives of at x = and x = 1 have opposite 
slopes, 

/(!)'(()) = (9) 

Therefore the perturbations f a under consideration have a singularity of order a in 
their derivatives at the end points of the unit interval, which are identified by the 
modulo one operation. 

Possible examples of such perturbations are polynomials or trigonometric func- 
tions. Such examples are 

f a (x) = x + ax(l -x) + 0{a 2 ) , (10) 
Ux) = x- -sin(7rx) +0(a 2 ) . (11) 

7T 



3. Hyperbolic regime. Although tp a ^ is a singular map, it has a unique abso- 
lutely continuous invariant measure, provided it is everywhere expanding, that is 



V'a{x) > 1. (12) 

The density of this measure has very interesting properties. Indeed, a direct con- 
sequence of the singularity of f a at the origin is that the iteration of a density of 
points initially uniform on the interval will yield a discontinuity which, provided a is 
small, is of the order of a. At the next iteration this discontinuity is rotated around 
the circle by £. Further iterations yield more and more discontinuities which, as 
long as C is irrational, will cover densely the whole interval. 
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Figure 1. Illustration of the map ip a £ denned by Eq. (0, with 
f a (x) the full form of Eq. i|TT)l given later in Eq. (|25[l. Here a = 1.5 
and C = 0.75. 



The invariant density is a fixed point of the Frobenius-Perron operator, 
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where we dropped the parametric dependences to simplify the notation. 
Let us consider the first order term in the expansion of p in a, namely 



p(x) = 1 + ap^\x) + 0(a 2 ) , 
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where we have taken the uniform density, p(x) = 1 to be the absolutely continuous 
solution of Eq. 1)13(1 when a — 0. An expression for follows from Eq. I|13(l : 



A 



< x < C, 



(15) 



i-C 



c < x < l. 



Let us comment here that p^ is generally not expected to be a well-behaved 
function of its arguments. This is readily seen if one notices that the two RHS of 
Eq. I|15|) evaluate to different values at x = £. This is due to the property Eq. @ 
that the derivative of is discontinuous at the origin. Moreover the discontinuity 
is of order of 5 = 0[f^ (0)]. Similarly, at x = 3£ mod 1 - the image of x = ( 
under <^o,Cj Eq. - another discontinuity appears, this time of order 5/2. In 
general, a discontinuity of order 5/2 n occurs at the nth iterate of £ under fox, i- e. 
a,t x — (J27=o m °d 1- For rational values of £, its binary expansion is a finite 
series : there exists a finite integer N such that £ = 5Zj=o ^i/2* +1 , w$ 6 {0, 1} and 
Wjv = 1. In that case, after N — 1 iterations of </?o,Ci C is mapped back on one of 
the iterates already visited, namely the one corresponding to the next largest index 
i for which oj,; = 1. Therefore, for rational £'s, the set of discontinuities is finite. 
However this is not so for irrational £'s. I n that case, the set of discontinuities is 
dense, even though of exponentially small amplitudes. 

Hence p^ 1 ' is a distribution, the formal expression of which can be be found by 
representing it as a Fourier series : 



k=l 



[A k (0 cos(27rfcx) + Sfc(C) sin(27rfea:)] 



(16) 



The coefficients A k (Q and B k (Q satisfy the recursion relations 

MO = ~9k sin(2^fcC) + A 2k (C) cos(27rfcC) - B 2k (C) sin(2^fcC), 

B k (0 = 9k cos(27rfcC) + A 2k (C) sin(2^fcC) + B 2k (() sin(2^fcC), 



with solutions 



Ak(0 = -^.92" fc sin[2(2™+ 1 -l) 7 rfcC] 

n=0 



where 



Bk(C) = ^52" fc cos[2(2" +1 -l)^C], 
g k = 8nk I dxf^ (^—^j cos(27rfc:r). 



(17) 
(18) 

(19) 
(20) 

(21) 
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It is not difficult to show that gu decays like 1/fe for k 3> 1. Indeed provided all 
the derivatives of exist at x = 0, the integral on the RHS of Eq. 1)21(1 can be 
written in the form of a series in l/k 2n (n integer) with bounded coefficients. 

Hence the series in Eqs. (I19I2U|) arc uniformly convergent trigonometric series, so 
that Ak and Bk are both continuous functions of £. We notice a striking similarity 
between the expressions of the Fourier coefficients Ak and Bk and Weierstrass 
functions |21j . which are continuous, nowhere differentiable functions. As stated 
in Appendix El below, the fractal dimension of these functions depends on the 
scaling properties of the coefficients in the trigonometric series. In particular, the 
dimension of the graphs of Ak and Bk defined by Eqs. I|19l20fl is unity, since the 
coefficients in the series decay like l/2 n . Figure (J2J) displays the first few Fourier 
modes Ak and Bk as functions of £, for k = 1, 2, 3. Note that modes of different fc's 
have similar shapes but decaying amplitudes as k increases. 

Consider again Eq. (|16|) . which, with the help of Eqs. 119I20|I . can alternatively 
be written in the form 

oo oo 

pW (x) = - J2 E sin {2nk [(2» +1 - l) C - x] } . (22) 

k=l n=0 

As we pointed out already, the dependence of p^ 1 ' upon x and C, is more complicated 
than the dependence of Ak,Bk on (. The set of its discontinuities is dense with 
respect to both arguments, even though p^ takes finite values everywhere. Nev- 
ertheless, the explicit formal functional form allows us to calculate some quantities 
of interest which involve p^ only under the integration sign. Such quantities are 
related to the Lyapunov exponent of the map and the transport properties of the 
same map "lifted" onto the real line. This is the subject of the next section. 

4. Lyapunov exponents and transport properties. The positive Lyapunov 
exponent can be computed up to second order in a with the help of Eqs. (|14I16|1 : 



A+(o,C) - 




With the expression of Bk(() given by Eq. I)2()[l. the summation in the last line 
can easily be computed numerically. 

Figure © shows the second order correction to A+ as a function of £. The fact 
that it is a symmetric function with respect to C ~ > 1 — C is a consequence of a 
symmetry of the map under ( — > 1 — £, a —> —a, x — > 1 — x and the fact that the 
positive Lyapunov exponent is a function only of the even powers of a. 

The stationary drift velocity measures the exchange of particles between the unit 
cells of the periodically extended version of f a x- ^ n this extension, based on the 
multi-baker map, points are sent to corresponding points in the unit intervals to 
their right and left after each iteration of the map. Points in the interval (/ a ((l — 
£)/2),/ a ((2 — C)/2)), mod(l), are sent one unit interval to the right, and assigned 
a velocity v = +1, while the remaining points in the interval (0, 1) are sent one 
unit to the left and assigned a velocity v = — 1, with identical behavior in each 
unit interval of the real line. This is, in fact, the projection onto the x-axis of the 
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FIGURE 2. From top to bottom and left to right : Ak(Q for k = 1, 
2 and 3, and Bk(C) tor k — 1, 2 and 3, numerically computed from 
Eqs. (|19I2U|) . where was taken from Eq. JTTJl 
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Figure 3. Absolute value of the second order correction to 
A+(a,£), Eq. Ij24(l . Here and in the next figure, g is derived from 
Eq. (f[J. 



nonlinear multi-baker map described in Appendix B (see Eq. H350 ). The expression 
for the stationary drift is now straightforward and yields 

r/-a((i-0/2) 
v(a, C) = dxp(x) 



/_ o ((2-0/2) /-l 

dxp(x) — / dxp(x). 

/_ o ((l-0/2) y/_ a ((2-C)/2) 

= ^ cos — -sin — 

oc 

+ EEt cos t( 2 " +2 " : ) ^ } + 0(fl3) ' (24) 

fc odd n=0 

Its graph is displayed in Fig. The drift velocity is an odd function of £ — 1/2. 
An interesting feature of Fig. J3J is that the drift velocity is negative in a small 
interval about £ = 1/2. This is another example of the negative current phenomenon 
discovered by Groeneveld in one dimensional maps, and discussed by Groeneveld, 
Klages, and coworkers ^JEIEI], an0 - ^ loosely corresponds to a negative electrical 
conductivity, where the current is in a direction opposite to that of the electric field, 
taken to be in the direction of C — I/ 2 - 
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FIGURE 4. Drift velocity to first order in a, Eq. (J2U- 

In Appendix IbI we present a generalization of our one-dimensional maps to two- 
dimensional maps, known as multi-baker maps. These maps have both a positive 
and a negative Lyapunov exponent. We believe that their sum obeys an identity 
relating the phase space contraction rate to the drift velocity, and is given as Eq. 
We do not yet have a rigorous proof of this identity, however. 

5. Non hyperbolic regime. Next, we describe some of the interesting features 
of the map ip a ,c f° r higher values of the field parameter. In Appendix[21 we present 
a classification of the non-linear maps f a satisfying the conditions listed in Sec. |2l 
We will restrict our attention to f a given by 



2 



7T 



TTX 



f a (x) = — arctan tan ( -77- ) e a , (25) 



2 7 



which is derived in Appendix B. As was mentioned in for £ = 0, the map 
loses its hyperbolicity for a > a c = log(2) and the origin becomes an attractive 
fixed point of the reduced map, so that, on a periodic lattice, the stationary state 
motion of particles is ballistic and moves toward decreasing n's. 

Things change quite dramatically as £ is increased. The origin is now mapped to 
£ and, for £ large enough, the orbit of the origin can eventually become hyperbolic, 
even though the slope of ip a ^ at x = would be less than one. The fixed point is 
not the origin x = anymore and it becomes attractive for larger and larger values 
of a as one increases £, i. e. a c (£), the value of the field parameter at which the 
fixed point becomes attractive, is a monotonically increasing function of £. 
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Moreover, for fixed £, stable periodic orbits can be observed in the range of a 
between log(2) and a c , as shown in Fig. JSJ. A property of these windows is that 
they all have as a fixed point for the value of a right before they bifurcate to 
chaotic windows. Based on this observation, it is possible to work out a condition 
for the presence of a window of stable period n trajectories, namely that 

<<(0) = 0, 

d 



dx 



<pW 



< 1, (26) 



x=0 

where the superscript n denotes the nth iterate. The system Eq. 112611 can be studied 
in order to find the minimal value of £ for which there exists such a stable period 
n orbit and the corresponding value of a at which it occurs. 

We illustrate this procedure for the simplest case of n = 2. From the top line of 
Eq. we find 

2tan(7rC/4) 



tan(7rC/4)] 



2 ■ 



(27) 



while the bottom line gives 



e > !| cos ^ . (28) 

~ 1 + cos(ttC) 

The minimal value of C satisfying Eqs. $Z7\ US} is C — 0-522873, for which a ~ 
1.005053. The existence of this periodic orbit is confirmed numerically. 

Higher n's are found to start appearing at lower values of C, which means that, for 
a fixed C, all the periods n greater than the first existing one will appear and they do 
so as a increases until it reaches its critical value a c . Such a succession of stable pe- 
riodic windows with increasing periods is referred to as a period adding-bifurcation, 
|15|. Here stable periodic windows alternate with chaotic ones. The similarity of 
the bifurcations observed here with period-adding bifurcations of reversible circle 
maps can be explained qualitatively if one notices that the two branches of (p a x 
have very different widths. Indeed the first critical point x = / a (l/2) gets closer 
to x = 1 as a increases, so that the first branch is close to a reversible circle map, 
whereas the second branch has a very steep slope. In this sense, for a fixed large 
a, as we vary £, we expect to see periodic orbits that remain on the first branch 
and thus correspond to stable windows, while the second branch accounts for the 
chaotic windows separating stable periodic windows. 

We illustrate this in Figs. (|5I7|I for the case of £ = 0.3, 0.5, and 0.7. In the 
first example, £ — 0.3, stable periodic windows are barely visible but one can check 
the existence of a period 6 window at a ~ 1.4449 whose width is O(10 -4 ). As C is 
increased to 0.5, windows get larger and, at a = 1.6492 a period 3 window starts 
that ends at a — 1.6881, then a period 4 window arises at a = 1.9678 that ends 
at a = 1.9985, and so on. But this is not the end of the story : besides the main 
sequence, there are a number of subsequences, such as In + 1, 3n + 1, 3n + 2, etc., 
some of which are clearly distinguishable on the third example, C = 0.7. The higher 
£, the greater the number of such subsequences. 

6. Conclusions and Discussion. In this paper we considered a class of two- 
parameter singular circle maps. One of these two parameters, denoted a, measures 
the intensity of the non-linearity. In the generalization of the non-linear baker map 
considered in Appendix iBl it plays the role of a thermostatted electric field. The 
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Figure 5. The bifurcation diagrams of <p a ^(x) for ( = 0.3 as a 
function of a. 



other parameter, denoted C, represents an external field, similar in some ways to a 
magnetic field. We found that the map exhibits both hyperbolic and non-hyperbolic 
behaviors, and that for small non-linearities, we could determine analytically the 
SRB measure of the fractal attractor to first order in a and the first non zero 
correction, of order a 2 , to the positive Lyapunov exponent. For large enough a we 
observed a transition to non-hyperbolic behavior with a period-adding bifurcation 
sequence leading to a stable fixed point for large enough a. 

Thus, the rotation parameter leads to a much more complicated dynamical be- 
havior of the map, including a region of "negative currents" similar to those found 
before by Groeneveld, Klages, and coworkers [1211141 ITU] , 

The results of this paper lead us to the following observation : The dynamical 
systems studied here are, in their hyperbolic regions, not Anosov diffeomorphisms. 
We find that the SRB measure of the fractal attractor is not a differentiable function 
of the parameter £, as it would be if the map were an Anosov diffeomorphism 
(i. e. did not have a singularity) \7\. Instead, it is a singular function of 
that parameter. As hyperbolic systems which are not Anosov occur frequently in 
models of physical systems it is not unreasonable to suppose that there will exist 
situations where non-differentiable SRB measures occur in physical systems. Our 
work suggests, at least, that this might very well be the case in Lorentz gases or 
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a 



Figure 6. Same as in Fig. JSJ for ( = 0.5. 

other discontinuous systems in the presence of external fields, where singularities 
of the flow would be rotated as the field is varied. 
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Appendix A. Weierstrass function. For A > 1 and 1 < s < 2, the Weierstrass 
function |2J can be defined as 

oc 

Wx, s : [0,1] — >R, t — > ^ A^- 2)fc sin(A fc i) . (29) 

k=l 

It is a nowhere differentiable continuous function |13j . 

As stated by Falconer [S] , provided A is large enough, the box counting dimension 
of the graph of this function is 

dim B Q{W x , s ) = s . (30) 
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Figure 7. Same as in Fig. JSJ for £ = 0.7. 



In general, the function 



F x , s (t) = Y\( s -V k g(\ k t) (31) 



k=l 



has box counting dimension s. This constitutes an upper bound for the Hausdorff 
dimension. However, it is commonly believed 3 that the Hausdorff dimension is 
exactly equal to s. 

The case s = 1, which is relevant to the Fourier modes given by Eqs. I|19I20[1 . 
is borderline. In this case, despite the nowhere differentiability of the Weierstrass 
function, it is understood that its Hausdorff dimension is equal to 1 pfl 117). 

Appendix B. Non-linear baker maps. Motivated by the analytical study of 
non-equilibrium stationary states of thermostated field driven systems, Gilbert et 
al. introduced in |llj the one parameter family of maps 

M a (n, x, y) = (1, 1, /„) o M ° (1, /„, 1) (n, x, y), (32) 

where 1 denotes the identity operator, acting cither on integers or on points of the 
unit circle, Mq is the usual multi-baker map |51 121)1 !§]. 

f (n-l,2x,%), < x < i 

M (n,x,y) = \ ,y + U ~x ^ *' ( 33 ) 

[ (ra + 1, 2x - 1, ), 5 < x < 1, 
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with variables n G Z and (x, y) £ [0, l] 2 . Here f a is a non-linear perturbation of the 
unit interval onto itself, which models the action of a thermostated external field, 



f a {x) — — arctan 

7T 



tan 



"KX 



(34) 



This model is a caricature of the two-dimensional periodic thermostated Lorentz 
gas with an applied electric field JHj Here the variable nx plays the role of the 
angle that the velocity makes with the direction of the electric field. 

We generalize Eq. (|32[) by introducing, in the multi-baker map, Eq. (|33|) . a new 
parameter, £, defined modulo 1, acting on the unit square, which shifts points by an 
angle Q in the rr-direction and acts on the y-direction so as to guarantee a weakened 
form of reversibility, see Eq. (|36[) below. 

The new multi-baker map, Fig. iJHJ, constructed so as to reduce to the multi- 
baker map, Eq. (|33[1 . for £ = and 1, reads : 



M QX (n,x,y) 



l,2x + C 



y 



2-C' 
"2 



< x < 



l,2x + C,^), 



i-C 

< y < C, 
i-C 



< a; < 



l,2x + C-l 
l,2x + C-2 
1,2a; + C- 2 



C<y< i, 
-C <T< 2-c 

0<y< 1, 

2^ < X < 1, 

< y < C, 

2-C 



(35) 



~2~ 



< x < 1, 



C<y< i. 
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Figure 8. M ox for C = 0.25. 
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The action of the time reversal operator T, T(x, y) = (1 — y, 1 — x), on this map 
can be seen to give 

T o M 0X oT = M-J_ ( . (36) 

Thus, under time reversal, the map (|35|) becomes the inverse of Mo : i_£ instead 
of the inverse of Mq^. This property is reminiscent of the effect of time reversal on 
an external magnetic field, which changes the sign of the field. 

We now define the two parameter map M a ^(n, x, y), the same way as M a in Eq. 
C2), by 

Max(n,x,y) = (1, 1, /„) o M £ o (1, f a , 1) {n,x,y), (37) 

Despite the lack of a precise reversibility, it is interesting to ask whether the the 
usual result relating the irreversible entropy production to the sum of the Lyapunov 
exponents still holds : 

5 = -(A+ + A_). (38) 

We believe it is the case here. This view is motivated by the symmetries present. 
However a definite answer based upon rigorous arguments has yet to be reached. 



Appendix C. Classification of time-reversible perturbations . In this ap- 
pendix we present a classification of the perturbations f a (x) , with fixed points at 
x = and x = 1, that satisfy the time-reversibility properties, Eqs. I|3l5f) . In order 
to find explicit forms of perturbations, it is natural to assume factorizability of the 
x and a dependences. Thus consider the two classes 

fa(x) = g- 1 o \g(x)q(a)] , (39) 
f a (x) = h^olh^+ria)} . (40) 

Here the parameter functions q and r are any functions of a with the respective 
properties, 

q(0) = 1, (41) 

q{a) n = q{na), (42) 

q(a)q(-a) = 1, (43) 

r(0) = , (44) 

nr(a) = r{na) , (45) 

r(a)+r(-a) = 0, (46) 

where n is a positive integer. In particular these conditions imply the property Eq. 
(@J for f a . We note that any function r of the second class defines a function q of 
the first class by exponentiation, q(a) = exp[r(a)]. The only function r satisfying 
Eqs. H44I46|I is a linear function r(a) = ca, where c is an arbitrary constant that can 
be taken to be c = — 1 without loss of generality, since a is allowed to take any real 
value. Likewise q(a) = exp(— a) is the only function that needs to be considered. 
The reversal symmetry, Eq. (JSJ) , is satisfied provided 

1 - ° [9(x)q(a)} - g- 1 o [g(l - x)/q(a)} , (47) 
1 - h- 1 o [h(x) + r{a)} = h- 1 o [h(l - x) - r(a)] . (48) 
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A class of functions verifying Eq. (|47|) is the set of functions g with the properties 

9~ l {x) = i-jT 1 ^), (50) 

Notice that the requirement that x = 0,1 be fixed points of f a implies that g is cither 
zero or plus or minus infinity at these points. Such examples are g(x) = tan(7ra/2), 
cot(7rcc/2), which define identical / a 's up to a change of sign of a. 
Similarly a class of functions verifying Eq. (|48|) is defined by 

h(x) = -h(l-x) , (51) 

h~ l {x) = l-h-\-x), (52) 

examples of which are h(x) = cot(7rx), cos(7ra;). Note however that cos(7rx) is not 
consistent with the requirement that x = 0, 1 remain fixed points of f a . Moreover 
the choice cot(7rx) yields a perturbation /„ with smooth derivative at the origin, 
/„(0) = /„(1) = 1. This can be easily seen to yield a differentiable density . 

We therefore argue that the choice of f a given by Eq. (|25|l is generic. Other 
forms of perturbations such as polynomial in powers of a can be constructed, but 
they do not have a complete form and can only be considered for values of a small 
enough with respect to the degree of the polynomial approximation. Thus the non- 
hyperbolic regime, where a becomes large, can only be studied with the example 
given in Sec. 
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